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In 2005, human coronavirus HKU1 (HCoV-HKUl) 
was isolated and identified from a 71-year-old 
man with pneumonia in Hong Kong. To identify 
and classify genotypes of HCoV-HKUl in Korea, a 
sensitive, specific, and quantitative real-time po¬ 
lymerase chain reaction (PCR) assay was devel¬ 
oped and analyzed the sequences of HCoV-HKUl 
isolated in Korea. A total of 1,985 respiratory 
specimens taken from patients with acute respi¬ 
ratory illness were tested for HCoV-HKUl from 
January 2007 to May 2008. The major clinical 
symptoms associated with HCoV-HKUl infection 
were examined statistically and sequence varia¬ 
tions of the RNA-dependent RNA polymerase 
(RdRp), spike, and nucleocapsid genes were also 
analyzed. Fifty cases (2.5%) HCoV-HKUl were 
identified by real-time PCR and viral loads ranged 
from 6.7 x 10 4 to 1.6 x 10 9 copies/ml. The clini¬ 
cal symptoms of HCoV-HKUl infection included 
rhinorrhea (72%), cough (64%), nasal congestion 
(56%), fever (32%), sputum (30%), sore throat 
(18%), chills (16%), postnasal discharge (14%), 
and tonsillar hypertrophy (10%). There was a sea¬ 
sonal distribution of HCoV-HKUl infection, peak¬ 
ing in winter and spring. Both genotypes A and 
B were detected but no recombination between 
them was found. This is the first report on the 
identification and genotyping of HCoV-HKUl as a 
causative agent of acute respiratory illness in 
Korea. The data suggest that at least two geno¬ 
types, A and B, of HCoV-HKUl with scattered 
silent mutations were circulating in Korea from 
2007 to 2008. J. Med. Virol. 85:309-314, 
2013. © 2012 Wiley Periodicals, Inc. 
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INTRODUCTION 

Human coronavirus (HCoV) 229E and HCoV-OC43 
account for 5-30% of cases of acute respiratory tract 
infection [McIntosh et al., 1969]. HCoV-NL63 was dis¬ 
covered in 2004 as a minor pathogen of HCoV infec¬ 
tions in humans [van der Hoek et al., 2004]. A new 
human coronavirus-HKUl (HCoV-HKUl) has been 
reported from a patient with pneumonia in Hong 
Kong. This species was classified into family Corona- 
viridae, genus Betacoronavirus, and further as a line¬ 
age A [Woo et al., 2005c]. 

The genome contains the genes for replicase poly¬ 
proteins la/lb (ppla/lb), hemagglutinin esterase 
(HE), spike (S), envelope (E), membrane (M), and nu¬ 
cleocapsid (N) in the same order as in other members 
of the genus Betacoronavirus [Pyre et al., 2007; Woo 
et al., 2009]. HCoV is associated frequently with com¬ 
munity acquired upper respiratory tract infections 
showing typical clinical symptoms including a runny 
nose, fever, coughing, and wheezing [Labret et al., 
2005; Esper et al., 2006; Gose et al., 2011J. However, 
other disease manifestations including bronchiolitis 
and pneumonia can occur [Bosis et al., 2007; Gerna 
et al., 2007; Vabret et al., 2007]. The respiratory 
symptoms of HCoV-HKUl infection are usually simi¬ 
lar to infections with other HCoVs, but HCoV-HKUl 
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can cause death in elderly people [Ren et al., 2011] 
and upper and lower respiratory tract infections in 
children [Woo et al., 2005b; Jin et al., 2012]. 

In this study, a Taq-Man ®-based real-time polymer¬ 
ase chain reaction (PCR) method that targets the 
HCoV-HKUl open reading frame (ORF) la and ORF 
lb genes with high sensitivity and specificity was 
developed and evaluated. Sequence information for 
the RdRp, S, and N genes of HCoV-HKUl was also 
analyzed to classify the genotype prevalence of the 
virus in Korea [Woo et al., 2005a; Woo et al., 2006]. 
This study is the first to identify and classify the 
genotypic predisposition of HCoV-HKUl as a causa¬ 
tive agent of acute respiratory illness in Korea. 

MATERIALS AND METHODS 
Clinical Specimens 

Throat swab specimens were collected from patients 
with an acute respiratory illness. Enrollment of 1,985 
patients was performed under the acute respiratory 
tract infection surveillance (ARI-Net) program in the 
Korea Centers for Disease Control and Prevention 
(KCDC) from January 2007 to May 2008. All clinical 
specimens were collected under regulation of Korea 
National Institute of Health institutional review 
boards (IRBs) and written informed consents were 
obtained from patients or their parent/guardian. 

Real-Time PCR Assay for HCoV-HKUl 

To design primer and probes for real-time PCR, full 
genome sequence of HCoV-HKUl was retrieved from 
GenBank accession no. NC006577. HCoV-HKUl ORF 
la and ORF lb gene were analyzed using Primer Ex¬ 
press software version 2.0 (Applied Biosystems, Foster 
City, CA) to allocate target sites which could be used 
as primers and probes. Total RNAs from 140 pi of 
clinical specimens were extracted using QIAamp viral 
RNA mini kit (Qiagen, Hilden, Germany) and 
directly used as templates. cDNA synthesis was per¬ 
formed using random primers and murine mammal 
leukemia virus (M-MLV) reverse transcriptase (Prom- 
ega, Madison, WI) in accordance with manufacturer’s 
instruction. Real-time PCR assays targeting the ORF 
la and ORF lb sequences were performed using 10 pi 
aliquots of mixtures containing 1 pi cDNA or a stan¬ 
dard plasmid, 5 pi Taq-Man l! universal PCR Master 


Mix (Applied Biosystems) containing ROX as a pas¬ 
sive reference dye, 900 nM of forward and reverse pri¬ 
mers and 250 nM of probe (Table I). Amplification 
and detection were performed using an ABI PRISM 
7900HT sequence detection system (Applied Biosys¬ 
tems) under the following conditions: uracilJV-glyco- 
sylase was activated at 50°C for 2 min, Taq 
polymerase was activated at 95°C for 10 min and 40 
cycles of amplification (95°C for 15 sec, 60°C for 

1 min) were run [Choi et al., 2008; Lu et al., 2008]. 
Positive values were decided when the copy number 
exceeded 10 copies/pl and the two target genes (ORF 
la and ORF lb) were detected simultaneously. The 
viral load was calculated and determined from mean 
threshold cycle (Ct) values of specimens and a stan¬ 
dard curve prepared from a known concentration of 
the DNA encoding the RdRp gene. All quantification 
experiments were performed using triple repetitions. 

Statistical Analysis of Clinical Symptoms 

The characteristics of clinical symptoms of HCoV- 
HKUl infection were investigated retrospectively 
from list of clinical records which were compiled up 
through ARI-Net. Regression analysis was performed 
with 95% confidential level to identify relationship be¬ 
tween HCoV-HKUl infection and clinical symptoms 
using the SAS program (version 9.2). 

Reverse Transcription (RT)-PCR and Sequencing 
of RdRp, S, and N Genes of HCoV-HKUl and 
Phylogenetic Analysis 

The RdRp, S, and N genes of HCoV-HKUl were am¬ 
plified by RT-PCR with specific primers (Table II). 
cDNA synthesis was performed as described above 
using M-MLV reverse transcriptase (Promega). The 
PCR assays were performed in a 50 pi reaction vol¬ 
umes containing 5 pi cDNA, 5 pi 10 x PCRbuffer, 

2 pi dNTPmix (final concentration of 400 mM of each 
dNTP), 1 pi SP Taq (Cosmo Gentech, Seoul, Korea), 
1 pi of each primer (10 pM), and nuclease-free water 
to 50 pi. The PCR reaction was carried out at 94°C for 

3 min, 35 cycles of amplification (30 sec at 94°C; 
30 sec at 48°C for RdRp, 58°C for S, 56°C for N; 30 sec 
at 72°C), and a final extension step at 72°C for 
10 min. The PCR products were run on a 2% agarose 
gel, stained with SYBR safe DNA gel stain dye (Invi- 
trogen, Carlsbad, CA) and visualized under UV light. 


TABLE I. Primer and Probe Sequences of HCoV-HKUl Used in Real-Time PCR Assay 


Gene 

Primer/Probe 

Sequence (5'-3') 

Position 

ORF la 

Forward 

GTTGGTTGTATGATGCGTTTGTTCT 

7,922-7,946 


Reverse 

TCTACAAATAAACTAGCATCAACATCATCGT 

7,971-8,001 


Probe 

FAM-CACGCTGTCCATCTCT-NFQ 

7,952-7,967 

ORF lb 

Forward 

CCTAACTGTGATGTGAGTGATGTCA 

16,516-16,540 


Reverse 

ACAAACCAAAGACCATACCATTCATAACT 

16,608-16,636 


Probe 

FAM-CATACCGCCCAAATAT-NFQ 

16,548-16,563 


FAM, 6-carboxylfluorescein; NFQ, non-fluorescent quenchers. 
J. Med. Virol. DOI 10.1002/jmv 
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TABLE II. Primer Sequences of HCoV-HKUl for Amplification of RdRp, S, and N Gene 

Primer 

Sequence (5'-3') 

Position 

RdRp-F 

RdRp-R 

S-F 

S-R 

N-F 

N-R 

GGG TAT GAA GTA TCA TCC TA 

GAT AAT CCC AAC CCA TAA GAA C 

AAC RYG GTG TTA TTA CTA 

AGA WGA TTG CAR AAA RCC AGA ACT 

ACT TGA RCG AAA YYA YCA AAC 

CGY AAA CCT AGT AGG GAT AGC TT 

14,433-14,452 

15,400-15,421 

23,756-23,773 

24,154-24,177 

28,409-28,429 

28,825-28,847 


Sequencing reactions were performed using BigDye" 
Terminators version 3.1 Cycle Sequencing Kits (Ap¬ 
plied Biosystems) using the manufacturer’s instruc¬ 
tions on a DNA analyzer (model 3730, Applied 
Biosystems). Phylogenetic tree construction was per¬ 
formed using the MEGA4 program (http://www.mega- 
software.net/) using the neighbor-joining method 
with Clustal W, and bootstrap values calculated from 
500 trees. 


Genbank Nucleotide Sequence 
Accession Numbers 

The following nucleotide accession numbers were 
used (http://www.ncbi.nlm.nih.gov/genbank/): refer¬ 
ence genes; AU884001, AY597011, DQ415899, 

DQ415911, DQ415912, DQ415914, DQ437616, 

EF507775, RdRp genes; JN234449, JN234450, 

JN234453, JN234454, JN234455, JN234456, 

JN234458, JN234459, JN234460, JN234461, 

JN234462, S genes; JN234465, JN234466, JN234468, 
JN234469, JN234470, JN234471, JN234473, 

JN234474, JN234475, JN234476, JN234477, N genes; 
JN234478, JN234479, JN234482, JN234483, 

JN234484, JN234487, JN234489, JN234490, 

JN234492, JN234493, and JN234494. 

RESULTS 

Development of Real-Time PCR 
for Detecting HCoV-HKUl 

The sensitivity and specificity of this real-time PCR 
were evaluated. Tenfold serial dilutions of the stan¬ 
dard HCoV-HKUl plasmid DNA ranging from 1 to 
1 x 10 10 copies/jil were used as templates for compar¬ 
ing the sensitivity between real-time PCR and con¬ 
ventional PCR. The detection limit of real-time PCR 
assays was as low as 10 copies/pl. The sensitivity of 
the real-time PCR assays used in this study was 
1 x 10 5 times higher than conventional PCR (Fig. 1). 
Cross reactivity was also evaluated with cDNA of in¬ 
fluenza viruses A/B, parainfluenza viruses 1/2/3, re¬ 
spiratory syncytial virus, human rhinovirus and 
enterovirus, DNA of adenovirus, human bocavirus, 
and respiratory pathogenic bacteria ( Haemophilus 
influenzae, Legionella spp., Streptococcus pneumoniae, 
Mycoplasma pneumoniae, and Chlamydophila pneu¬ 
moniae). No non-specific amplification was detected in 
either assay system (data not shown). 


Detection and Clinical Characteristics 
of HCoV-HKUl Infection 

Fifty HCoV-HKUl-positive cases were detected 
(2.5%) out of 1,985 clinical specimens using real-time 
PCR assays targeting ORF la and ORF lb. The 
amount of RNA detected in positive specimens ranged 
from 6.7 x 10 4 to 1.6 x 10 9 copies/ml. Based on analy¬ 
sis of multiple alternative clinical records, rhinorrhea 
(72% of patients), cough (64%), and nasal congestion 
(56%) were the major symptoms of patients infected 
with HCoV-HKUl. The other basic respiratory symp¬ 
toms were fever (32%), sputum discharge (30%), sore 
throat (18%), chills (16%), postnasal discharge (14%), 
and tonsillar hypertrophy (10%; Table III). There 
were no positive cases associated with severe respira¬ 
tory symptoms such as pneumonia and bronchiolitis. 
Seasonal distribution of HCoV-HKUl infection was 
found, with the main prevalence from early winter to 
late spring with an 8% mean positive rate (Fig. 2). 

Sequencing and Phylogenetic Analysis 
of HCoV-HKUl 

RdRp, S, and N genes and the putative amino acid 
sequences of 11 cases, of which the target genes were 
amplified by conventional PCR, were analyzed. Eight 
cases were HCoV-HKUl genotype B (73%) and three 
were HCoV-HKUl genotype A (27%) when queried 



Fig. 1. Comparison of the sensitivity of the real-time polymerase 
chain reaction (PCR) assay for the human coronavirus HKU1 
(HCoV-HKUl) ORF lb gene with conventional PCR. The RT-PCR 
assay detected up to 10 copies/|xl of standard DNA, whereas the de¬ 
tection limit of conventional PCR was 1 x 10 6 copies/pl. 
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TABLE III. Proportion of Clinical Symptoms of HCoV- 
HKU1 Infected or Non-Infected Patients 



No. (%) of 
HCoV-HKUl a 

No. (%) of 
non-HCoV-HKU l a 

Rhinorrhea 

36 (72.0) 

1,547 (79.9) 

Cough 

32 (64.0) 

1,318 (68.1) 

Stuffy nose 

28 (56.0) 

1,215 (62.8) 

Fever 

16 (32.0) 

622 (32.1) 

Sputum 

15 (30.0) 

493 (25.5) 

Sore throat 

9 (18.0) 

267 (13.8) 

Chill 

8 (16.0) 

230 (11.9) 

Postnasal discharge 

7 (14.0) 

332 (17.2) 

Tonsillar hypertrophy 

5 (10.0) 

114 (5.9) 

Throat flare 

4 (8.0) 

235 (12.1) 

Otalgia 

3 (6.0) 

19(1.0) 

Pain 

3 (6.0) 

52 (2.7) 

Hoarseness 

2 (4.0) 

78 (4.0) 

Eye discharge 

1 (2.0) 

37 (1.9) 

Otarrhea 

1 (2.0) 

10 (0.5) 

Wheeze 

0 (0.0) 

37 (1.9) 

Dyspnea 

0 (0.0) 

38 (2.0) 


“Clinical symptoms were scored as multiple-alternative manner. 


through a public database [Benson et al., 2005]. Phy¬ 
logenetic analysis also indicated that HCoV-HKUl 
isolates in this study consisted of genotypes A and B, 
whereas genotype C was not identified (Fig. 3). 

DISCUSSION 

Since HCoV-HKUl was first identified in a patient 
with pneumonia in Hong Kong, clinical manifesta¬ 
tions and epidemiological features of infection with 
this virus have been reported from several countries 
[Labret et al., 2005; Esper et al., 2006; Gose et al., 
2011]. Chung et al. [2007] also tried to identify the 
virus in patients with acute expiratory wheezing in 
Korea but evidence of HCoV-HKUl infection was not 
found. These negative results arose from several limi¬ 
tations such as an insensitive diagnostic strategy, se¬ 
lection of samples out of season as well as the 
relatively small number of patients studied. To cir¬ 
cumvent these limitations, a quantitative real-time 
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Fig. 2. Seasonal distribution of HCoV-HKUl infection from 
January 2007 to May 2008. The HCoV-HKUl prevailed from early 
winter to late spring in Korea. Bars indicate the monthly number of 
clinical specimens. The line represents the monthly proportion of 
positive cases as a whole. 
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Fig. 3. Phylogenetic analysis of the RdRp, S, and N genes of 
HCoV-HKUl. Eleven positive cases were classed as genotype A 
(n = 3) or B (n = 8). There was no HCoV-HKUl genotype C detected 
in Korea during the study period. Each shadowed number indicates 
the specimen identification number, and reference sequences are 
marked as GenBank accession numbers. Phylogenetic trees are 
shown of (A) the RdRp gene (897 neucleotide (nt)), (B) the S gene 
(309 nt), and (C) the N gene (425 nt). Primer sequences were exclud¬ 
ed from phylogenetic analysis. 


PCR was developed and targeted a study group of 
1,985 throat swab specimens collected from patients 
with an acute respiratory illness from January 2007 
to May 2008. 

The aim of the real-time PCR assays developed in 
this study was to detect both of the HCoV-HKUl ORF 
la and ORF lb genes. The assays provided enhanced 
sensitivity and specificity, which could permit quanti¬ 
tative analysis of HCoV-HKUl in clinical specimens. 
The ORF lb assay gave more dynamic Ct values than 
the ORF la assay. Therefore, the ORF lb assay could 
be more applicable for clinical diagnosis using respira¬ 
tory specimens from patients with acute respiratory 
tract infections caused by HCoV-HKUl. 

In this study, HCoV-HKUl was found in 2.5% of 
patients (50 cases) with acute respiratory illness sam¬ 
pled from January 2007 to May 2008 in Korea. A sim¬ 
ilar incidence was reported previously in Hong Kong 
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[Woo et al., 2009] and China [Jin et al., 2010] but an 
apparently higher rate was noted in Chinese study be¬ 
cause of different clinical practices. Unlikely, appar¬ 
ent lower incidence (0.3%, n = 4,181) of HCoV-HKUl 
in patients with acute respiratory tract in Hong Kong 
was reported by Lau et al. [2006]. Relatively high inci¬ 
dence of HCoV-HKUl in Korea might have resulted 
from difference of target study group and/or sensitivi¬ 
ty of detection strategy. Even though the quantified 
viral loads varied from 6.7 x 10 4 to 1.6 x 10 9 copies/ 
ml in the respiratory disease specimens used, there 
were no significant correlations between viral load 
and specific clinical features. Similarly, there were no 
statistically significant relationships between the 
HCoV-HKUl infection and any particular clinical 
symptoms. As for seasonal prevalence, the HCoV- 
HKU1 infections prevailed in the winter and spring in 
Korea and this result was consistent with the infec¬ 
tion rates of HCoV-OC43, 229E and with other reports 
[Lau et al., 2006; Domiquez et al., 2009; Gaunt et al., 
2010]. One limitation of this study was that because 
clinical specimens were chosen from negative cases of 
acute respiratory illness in the ARI-Net laboratory 
surveillance system [Chun et al., 2009], dual or multi¬ 
ple infections of HCoV-HKUl with other respiratory 
viruses could not be detected. 

To identify whether there was any recombinant 
HCoV-HKUl isolated from Korean patients, the nu¬ 
cleotide sequence information of the RdRp, S, and N 
genes was analyzed. Similar to previous study 
reported by Lau et al. [2006] in Hong Kong, our phylo¬ 
genetic analysis revealed that 11 cases of HCoV- 
HKUl isolated in Korea comprised two genotypes: A 
or B. When single reciprocal nucleotide exchanges 
were analyzed between the RdRp, S, and N genes, no 
recombination could be identified. Instead, several 
point mutations were observed at the nucleotide level, 
resulting in a transition (C to T) or a transversion (C 
to T or T to C) in the RdRp, S, and N genes. However, 
these were silent mutations so that no changes oc¬ 
curred at the amino acid level. Interestingly, a silent 
continuous mutation in the N gene was identified 
from 6/8 isolates of HCoV-HKUl genotype B. Thus, 
there might have been a unique strain of HCoV- 
HKUl prevailing in Korea at that time. 

This is the first report that HCoV-HKUl was asso¬ 
ciated with acute respiratory illness in Korea. Infec¬ 
tions with HCoV-HKUl peaked from winter through 
to spring in Korea during 2007-2008. Sequencing 
and phylogenetic analysis revealed that the 11 
HCoV-HKUl isolates belonged to two distinct geno¬ 
types A and B based on the RdRp, S, and N genes. 
Unlike a previous study, no gene recombination was 
identified in HCoV-HKUl. Instead, a conserved si¬ 
lent mutation was observed in the N gene of the pu¬ 
tative Korean strain of HCoV-HKUl genotype B 
from six specimens. Intensive molecular epidemiolog¬ 
ical investigations would be required to elucidate any 
possible changes in viral tropism caused by such 
mutations. 
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